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Abstract 

We consider the X-cell multiple-input multiple-output (MIMO) interfering multiple-access channel (IMAC) 
with time-invariant channel coefficients, where each cell consists of a base station (BS) with M antennas and N 
users having L antennas each. In this paper, we propose two opportunistic interference alignment (OIA) techniques 
utilizing multiple transmit antennas at each user: antenna selection-based OIA and singular value decomposition 
(SVD)-based OIA. Their performance is analyzed in terms of user scaling law required to achieve KS degrees- 
of-freedom (DoF), where S(< M) denotes the number of simultaneously transmitting users per cell. We assume 
that each selected user transmits a single data stream at each time-slot. It is shown that the antenna selection- 
based OIA does not fundamentally change the user scaling condition if L is fixed, compared with the single-input 
multiple-output (SIMO) IMAC case, which is given by SNR^~ 1 - ) ' S , where SNR denotes the signal-to-noise ratio. 
In addition, we show that the SVD-based OIA can greatly reduce the user scaling condition to SNR^~ 1 - ) ' S-i+1 
through optimizing a weight vector at each user. Simulation results validate the derived scaling laws of the proposed 
OIA techniques. The sum-rate performance of the proposed OIA techniques is compared with the conventional 
techniques in MIMO IMAC channels and it is shown that the proposed OIA techniques outperform the conventional 
techniques. 
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I. Introduction 

Interference management is a crucial problem in wireless communications. Over the past decade, there 
has been a great deal of research to characterize the asymptotic capacity inner-bounds of interference 
channels (ICs) using the simple notion of degrees-of-freedom (DoF), also known as multiplexing gain. 
Recently, interference alignment (IA.) [2]-[10] has emerged as a fundamental solution to achieve the 
optimal degrees-of-freedom (DoF£| in several IC models. The conventional IA technique for the K- 
user IC [2] and the i^-user X channel [9], [10] is based on several strict conditions as follows. Time, 
frequency, or space domain extension is required to render the channel model multi-dimensional. To 
this end, channel randomness, i.e., time-varying or frequency-selective channel coefficients, is needed. 
Moreover, an arbitrarily large size of the dimension extension is needed for K greater than 3, which 
results in an excessive bandwidth usage is required for the decoding of one signal block [6]. In addition, 
global channel state information (CSI) is needed at all nodes [2], [3], [8], [11], [12]. 

For the interfering multiple-access channel (IMAC) consisting of K cells, where each cell is composed 
of iV users and a single base station (BS), Suh and Tse developed a new IA scheme to characterize the 
DoF achievability of the i^-cell IMAC [6] allowing the rank of the interference space to be larger than 
one. The underlying idea of the IA is to align the interference to the desired interference spaces at the 
receivers by exploiting diversity (i.e., randomness) in any resource domain. The scheme proposed in [6] 
utilized the user domain resource for the IA in the IMAC. This IA scheme based on the user diversity leads 
to two interesting results. First, the DoF of the interference-free network, given by K, can be achieved 
as iV increases. Second, the size of the time/frequency domain extension is greatly reduced. Specifically, 
the finite size of the extension is given by O(N), which is sufficient to operate for given N. However, 
arbitrarily large N is needed to attain K DoF, which results in an infinite dimension extension in the end. 
Thus, time-varying or frequency- selective fading is still required for this scheme. 

Recently, the concept of opportunistic interference alignment (OIA) was introduced in [13]— [17], for 
the .fT-cell iV-user single-input multiple-output (SIMO) IMAC with time-invariant channel coefficients, 
where each base station (BS) has M antennas. In the OIA technique, opportunistic user scheduling 
is combined with the spatial domain IA to align the interference to predefined interference spaces at 
each BS by exploiting multiuser diversity. Although several studies independently addressed some of the 
aforementioned practical problems of the conventional IA technique [6], [18]-[20], the OIA technique 
resolved these practical issues simultaneously. The OIA scheme employs the spatial domain IA only with 
the aid of opportunistic user scheduling and thus operates with a single snapshot without any dimension 
extension. The purpose of the OIA-related work [13], [15], [16] is not only to maximize the DoF as in the 
conventional schemes, but also to characterize the trade-off between the achievable DoF and the number 
of users required. It was shown in [16] that the OIA scheme achieves KS DoF if iV scales faster than 
SNR (^-i)s in a nign SNR regime, where S{< M) is the number of selected users in each cell. 

In this paper, we introduce an OIA for the i^-cell MIMO IMAC with time-invariant channel coefficients, 
where each cell consists of one BS with M antennas and N users having L antennas each. Inheriting 
the basic OIA principle [13], the proposed OIA operates with local CSI at the transmitteiH, no inter-user 
or intercell coordination (i.e., distributed scheduling metric calculation), no dimension extension, and no 
iterative processing. In [21], the outer bound on the DoF of the MIMO IMAC with time-invariant channel 
coefficients was characterized, and necessary conditions for M and L needed to achieve the optimal DoF 
were derived with global CSI at all nodes. However, the main goal of the proposed OIA is to characterize 
a trade-off between the achievable DoF and the number of users required in the MIMO IMAC with 
arbitrary M and L. That is, the focus is on studying the user scaling law needed to achieve the target 
DoF, given by KS, which is optimal if S = M. Scaling conditions required to achieve target performance 
have a great impact in providing the convergence rate to the target performance with respect to considered 

'The optimal DoF denotes the maximum achievable DoF for given channel, which is proved by the converse proof. 
2 In interference channels, the local CSI at the transmitter denotes the information of the channels from the transmitter to all receivers, 
i.e., its own transmit links [5]. 
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system parameters, thus yielding an intuitive performance measure. For instance, it is common in MIMO 
systems to evaluate limited feedback schemes by analyzing the relationship between the codebook size 
scaling and the rate-loss [22], [23], and the concept has been applied also to MIMO ICs [24], [25]. 

In the downlink cellular IC, user scaling laws were developed for the OIA [21], [26] and for the 
opportunistic interference management with limited feedback [27], [28]. These schemes cannot be easily 
extended to the IMAC, because there exists a mismatch between generating interferences at each user 
and interferences suffered by each BS from multiple users, thus yielding the difficulty of user scheduling 
design. 

More specifically, we propose the following two types of OIA: antenna selection-based OIA and singular 
value decomposition (SVD)-based OIA. We then derive the scaling law for required N with respect to 
SNR, under which KS DoF can be achieved. In the proposed schemes, each selected user employs 
transmit beamforming to mitigate the leakage of interference (LIF) it generates. While the alignment was 
performed only through user scheduling in the SIMO case, the transmit beamforming is used for the 
MIMO OIA to perform the spatial domain IA along with opportunistic user scheduling. Moreover, the 
additional effort for the feedback of the weight vector from each selected user to the corresponding BS 
is in general required compared to the SIMO case, except for the proposed antenna selection-based OIA. 

We show that for the antenna selection-based OIA, where the best transmit antenna is selected at each 
user, required N scales as L^SNR^ -1 * 15 . Thus, the user scaling condition with respect to SNR does 
not fundamentally change, compared with the SIMO IMAC case [16], if L is a constant independent 
of N. However, the sum-rate gain of the antenna selection-based OIA over the SIMO OIA increases as 
L grows, whereas no additional feedback is required. For the SVD-based OIA, each user designs the 
weight vector that minimizes the leakage of interference (LIF) using SVD-based beamforming. We show 
that the SVD-based OIA can greatly reduce the user scaling condition to SNR < - X_1 ^ 5_L+1 with the help 
of the high-rate feedback. Our schemes are compared with the existing IA schemes for multiuser ICs, 
and computer simulations are provided to validate the derived scaling laws. From this study, besides the 
fundamental trade-off between the user scaling condition and the achievable DoF, we examine that in the 
MIMO IMAC, there also exists a trade-off between the amount of feedback for the weight vectors and 
the user scaling condition. 

The organization of this paper is as follows. Section QI] describes the system and channel models of 
MIMO IMAC. The proposed the MIMO OIA scheme is presented in Section Unl Both DoF achievability 
analyses and user scaling laws are provided in Section [IV] The proposed scheme is compared with the 
existing MIMO uplink schemes as well as the converse proof in Section [V] Section [VI] provides simulation 
results and Section IVIII concludes the paper. 

Notations: C indicates the field of complex numbers. (-) T and (-) H denote the transpose and the 
conjugate transpose, respectively. 

II. System and Channel Models 

Let us consider the time-division duplex (TDD) fT-cell MIMO IMAC, as depicted in Fig. [0 Each cell 
consists of a BS with M antennas and N users, each with L antennas. The number of users selected to 
transmit uplink signals in each cell is denoted by S < M. It is assumed that each selected user transmits 
a single spatial stream. To consider nontrivial cases, we assume that L < (K — 1)S + 1, because all the 
inter-cell interference can be completely canceled at the transmitters otherwise!!. The channel matrix from 
user j in the z-th cell to BS k (in the k-th cell) is denoted by B 1 ^ e C MxL , where i, k e fC = {1, . . . , K} 
and j 6 A/" = {1, ... , N}. Time-invariant frequency-flat fading is assumed, i.e., channel coefficients are 
constant during a transmission block, and channel reciprocity between uplink and downlink channels is 
assumed. From pilot signals sent from all the BSs, user j in the i-th cell can estimate the channels H^, 
k = 1, . . . , K, utilizing the channel reciprocity, i.e., the local CSI at the transmitter. Without loss of 

3 The case where L > (K — 1)S + 1 and where each selected user transmits multiple spatial streams is discussed at the end of Section 
IIV-BI (see Remark [3} and also in Section [V-BI with the comparison to the existing schemes. 
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generality, the indices of selected users in every cell are assumed to be (1, . 
defined by 



DoF 



lim 

SNR->oo 



Ef=iE=i^ [ 



logSNR 

where denotes the achievable rate for user j in the i-th cell. 



S). The total DoF are 



(1) 



III. Proposed OIA for MIMO IMAC 

We first describe the overall procedure of the proposed OIA scheme for MEMO IMAC, and then derive 
the achievable sum-rate and present the geometric interpretation of the proposed scheme. 



A. Overall Procedure 

1) Initialization (Reference Basis Broadcast): The interference space for the interference alignment 
at the k-th cell is denoted by Q fc = [q^i, . . . , qk,M-s]> where (\ k ,m G C Mxl is the orthonormal basis, 
k G /C, m = 1, . . . , M — S. BS k independently generates q fc m from the isotropic distribution over the 
M-dimensional unit sphere. For given Q fc , BS k also calculates the null space of Q^, defined by 



U fc = [u M , . . . , u k ,s] - null(Q A 



(2) 



where u kf i G C is the orthonormal basis, and broadcasts it to all users prior to the communication. The 
interference basis Q k can be chosen arbitrarily such that Q k is full rank. A simple way to maximize the 
performance of the ZF equalization at the BS, which will be discussed in the sequel, would be choosing 
M — S columns of the left or right singular matrix of any M x M matrix as Q k and choosing the rest 
of the S columns as XJ k . If S = M, then can be any orthogonal matrix. Note that the calculation and 
broadcast of is required only once prior to the communication as is determined only by M and 
S. 

2) Stage 1 (Weight Design and Scheduling Metric Feedback): Let us define the unit-norm weight 



vector at user j in the i-th cell by w l '' jj , i.e., w 



<•>.;] 



1. Two different methods to design w^ Ml shall be 



presented in Section JV] along with the corresponding user scaling law. From the notion of and H^, 
user j in the i-th cell calculates its LIF, which is received at BS k and not aligned at the interference 
space Qfc, from 



Proj 



UfHL IJJ w 



-LQ* 



H 



S3. 



W 



(3) 
(4) 



where i G K,, j G AT, and k G K, \ i = {1, 



+ K}. The scheduling metric of user j in the 



i-th cell, denoted by r^-^, is defined by the sum of LIFs, which are not aligned to the interference spaces 
at neighboring cells. That is, 



K 



7/ 



hi] 



E 

=l,k^i 



(5) 



All the users report their LIF metrics to corresponding BSs. 

3) Stage 2 (User Selection): Upon receiving users' scheduling metrics in the serving cell, each BS 
selects S users having smallest LIF metrics. Note again that we assume without loss of generality that 
user j, j — 1, . . . , S, in each cell have the smallest LIF metrics and thus are selected. Subsequently, user 
j in the i-th cell forwards the information on to BS i for coherent decoding. 
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4) Stage 3 (Uplink Communication): The transmit signal vector at user j in the i-th cell is given by 
W M X M 5 where is the transmit symbol with unit average power, and the received signal at BS i can 
be written as: 



s 



desired signal 
K S 

r[k,m] [k,m] Jfc.m] 



inter-cell interference 



where z, G C A/Xl denotes the additive noise, each element of which is independent and identically 
distributed complex Gaussian with zero mean and the variance of SNR -1 . As in SIMO IMAC [13], [16], 
the linear zero-forcing (ZF) detection is applied at the BSs to null inter-user interference for the home 
cell users' signals. From the notion of H 4 and w^', BS i obtains the sufficient statistics for parallel 
decoding 

r<= [r i>1 ,...,r i!S } T ±F i H Ufy i , (7) 

where Uj is multiplied to remove the inter-cell interference components that are aligned at the interference 
space of BS i, Qj, and Fj e C SxS is the ZF equalizer defined by 



[f»,i> • • • ) f;,s] 



U i ^HP 1 w^ 1 ],...,U i H Hf' 51 w^ 



(8) 



For a comprehensive overview, the overall sequential procedure is illustrated in Fig. |2] Note that we 
assume low-rate perfect information exchanges for Stages 1-3, such as feedback of the scheduling metric, 
broadcast of user selection information, feedforward of weight vector information, as in [16], [23]-[26], 
[29]. 



B. Sum-Rate Calculation 

From fT), the jth spatial stream, r^-, is written as 



+ fu B <, (9) 



where z\ = Uj Zj. Thus, R [l ' n is given by 

RM = log + SINR [IJ - 



Mi+ |, f )■ do) 

\\H,j || + 

where SINR'*'^ denotes the signal-to-noise-and-interference ratio of the user j in the i-th cell and Ijj is 
the sum-interference defined by 



k=l,k^i m=l 



SNR. (11) 
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C. Geometric Interpretation 

If S < M and the interference from user m in the A;-th cell to BS i is aligned to Qj, i.e., 

Hf^w^espantQ,], (12) 

then it is nulled in because U^Hf ' m W fc,m l = 0, i.e., ^| fc ' m ' = 0. If S = M, the LIF metric is simplified 



tO 7] l k J 



2 



In this case, no IA is conducted and only the opportunistic interference nulling 

(OIN) is performed as in the OIN for the SIMO IMAC [16]. We do not separately describe this OIN 
mode, as it can be taken into account by the OIA framework. 

Figure |3] illustrates the proposed MIMO OIA for K = 2, M = 3, and S — 2. The interference terms 
H2 w' 1,1 ' and H2 w' 1,2 ' should be aligned to the interference space c^i at BS 2, while we only require 
for the signal vectors H^w' 1,1 ! and H^'^w' 1 ' 2 ] to be distinguishable at BS 1. Similarly, Hf'^w' 2,1 ! and 
Hf' 2 W 2 ' 2 l should be aligned to q^i at BS 1, while H^'^w^ 2 ' 1 ' and Hjf^W 2 ' 2 ] need to be distinguishable 
at BS 2. The main task of the achievability proof is to show that r/[ 1 '^ can be made arbitrarily small for 
all cross-links through opportunistic scheduling and beamforming, which proves that the IA conditions 
(fl"2l hold true almost surely for all i G /C, m G S = {1, . . . , S}, and k G K \ i. Note that for given 
w' 1 '^, the signal vectors at each BS are distinguishable, since the channel coefficients are generated from 
continuous distributions. Therefore, in such case, the DoF of KS is achievable. 

IV. DoF Achievability 

In this section, we present two different beamforming strategies to design at each user, and 

characterize the DoF achievability for each strategy in terms of the user scaling law. 

A. Antenna Selection 

In the antenna selection-based OIA, only one transmit antenna is selected to transmit at each user, i.e., 
W M £ | ei) . . . ( ex}, where denotes the Z-th column of the (L x L) -dimensional identity matrix. Let us 
denote the Z-th column of by hj*'^, Z G {1, . . . , L}. Then, user j in the i-th cell chooses the optimal 
weight vector as wjjjj = e^^, where the index is obtained from 



K 

j) = argmm || u fe Hh 
Then, the corresponding scheduling metric is given by 

K 



h3i 
k,l 



(13) 



k,l(i,j) 



(14) 



and is reported to BS i. Since the l(i, ?')-th column of the channel matrix, h », . , is the effective channel 

i,l{i,3) 

vector at BS i, the feedback is not needed if user j in the i-th cell transmits the uplink pilot to BS i only 
through the j ) -th antenna after it is selected to transmit. 

The following theorem establishes the DoF achievability of the antenna selection-based OIA. 

Theorem 1 (User scaling law: Antenna selection-based OIA): The antenna selection-based OIA with 
the scheduling metric ([141 ) achieves 

DoF > KS (15) 

with high probability if 

N = u (l- 1 SNR ( - k - 1)s ) , (16) 
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where a function f(x) defined by f(x) = u(g(x)) implies that lim^^ = 0. 

Proof: See Appendix HI ■ 
Note that in the SIMO IMAC, the OIA scheme achieves the DoF of KS if N = co (SNR (i ^ 1)5 ) [16, 
Theorem 1]. Thus, the antenna selection-based OIA does not fundamentally change the user scaling if L 
is fixed. Note that however, the user scaling condition is reduced even without any additional feedback, 
compared to the SIMO case, if L scales with respect to SNR. The following remark discusses the 
cooperative feature the opportunistic gain obtained from the user and antenna diversity in the antenna 
selection-based OIA. 

Remark 1: If L scales faster than SNR^ L , where is a positive scalar, then the user scaling condition 
to achieve the DoF of KS is given by iV = u (sNR ( *~ 1)s-1 k ) . If i/> L = (K-1)S, then the DoF of KS is 
obtained with high probability for any iV > S. In such case, the opportunistic gain is sufficiently obtained 
only through the antenna diversity. In other words, the opportunistic gain can be achieved cooperatively 
from the user and antenna diversity. 

Now as a corollary to Theorem 1 in [16], we discuss the upper-bound on the user scaling law with 
the antenna selection by considering the general case where more than one transmit spatial stream are 
allowed at each user. 

Corollary 1: Suppose that user j in the i-th cell selects rv h ^ transmit antennas with smaller LIF metrics, 



where the /-th antenna's LIF metric is given by J2f=i k^i Ufc^hj^ J , / G {!,...,£,}. Then, the general 



H u {ij] 



2 



users with 



smaller sum-LIF metrics, achieves 



antenna selection-based OIA, in which BS i selects Si 
KS DoF with high probability if N = u (L~ 1 SWi iK '' 1)s \ and if n& j] and Si are chosen such that 

S = J2'jLi i € JC, and such that the sum-LIF of the selected S spatial channels is minimized at each 
cell. 

Proof: Since the considered scheme is equivalent to selecting S spatial channels (transmit antennas) 
with smaller LIF metrics amongst NL spatial channels, which is also equivalent to the SIMO OIA with 
NL users, the proof is immediate from [16, Theorem 1]. ■ 
Remark 2: In the general antenna selection approach, the optimal should be determined to find 
the best S spatial channels, which in general requires a joint optimization with global CSI or L times 
increased the feedback phases for each user to feed back all individual LIF metrics for L antennas. 
Surprisingly, Theorem [T] indicates that the antenna selection-based OIA with single spatial stream at each 
user is enough to achieve the same result, in which the scheduling metric is calculated at each user using 
local CSI without any cooperation or additional feedback. It is more surprising that the selection of the 
best one out of the L spatial channels at each user does not degrade the diversity gain in terms of the user 
scaling law, compared to the selection of the best S out of NL spatial channels. The result is encouraging, 
since we can expect the same benefit of increasing N to NL by simply increasing the number of antennas 
at the users. 



B. SVD-Based OIA 



In the SVD-based OIA, each user finds the optimal weight vector that minimizes its LIF metric. The 
same beamforming technique was also considered in [30], [31] for the MEMO IMAC, however, our focus 
is to derive a user scaling law and thereby to analytically examine the relationship between the number 
of users and the beamforming techniques used. 

The LIF metric for the SVD-based OIA is defined by 



VSVD 



K 

E 



(17) 
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where G^ G R -^ SxL is the stacked cross-link channel matrix, defined by 

T / „ r. . n \ T 



(u^Hp 1 ) , . . . , (u^HfcS 



Qlhi] A 

(U i+1 *H^) ,...,( U K »H K 

Let us denote the SVD of G [i ' j] as 



T / \ T 

H-rj[hj\\ In H-rjlhi] 



(18) 



G [*>•?] = f2'* ,J ''s'* ,5 ''v' i '^ H (19) 

where G C^" 1 ) 5 ^ and G C LxL consist of L orthonormal columns, and S [i ' i] = diag (a\' j \ ...,a 
where a\ > • ■ ■ > a I . Then, it is apparent that the optimal is determined as 

w^ ] D = arg min 1 1 G v 1 1 2 = vf ' j] , (20) 
where is the L-th column of V^l With this choice the LIF metric is simplified to 

rg& = 4' J]2 - (2D 



^svd va l ues among N users than the rest. To construct Fj defined in © at BS i for given selected user 
j, i — 1, . . . , K, j — 1, . . . , S, the information of w!^ needs to be known by BS i through the feedback 



All the users report their LIF metrics to the corresponding BSs and BS i selects S users with smaller 
77svd values among users 
j,i = l,...,K,j = l,...,l 
with a sufficiently high rate. 

At this point, we introduce a useful lemma for the polynomial expression of the CDF of 7?!^. 
Lemma 1 (CDF ofrj^): The CDF of r/|^, denoted by F a (x), can be written as 

F a (x) = ax^ s - L+l + o {x^ s - L+l ) , (22) 

for < x < 1, where a is a constant determined by K, S, and L. 

Proof: See Appendix HB ■ 

Now the following theorem establishes the DoF achievability of the SVD-based OIA. 

Theorem 2 (User scaling law: SVD-based OIA): The proposed SVD-based OIA scheme with the schedul- 
ing metric (121) achieves 

DoF > KS (23) 

with high probability if 

N = u (SNR (X " 1)5 - L+1 ) . (24) 
Proof: See Appendix [nil ■ 
Therefore, unlike the antenna selection, the SVD-based OIA fundamentally lowers the power of the SNR 
scaling condition required to achieve the DoF of KS. Note that however, this reduced scaling is achieved 
at the cost of the sufficiently high-rate feedback of Wgy^ from all the selected users to associated BSs. 
Noting that the antenna selection-based OIA needs no feedback, the antenna selection- and SVD-based 
OIA schemes are the two extremes of the trade-off between the feedback amount and the user scaling 
condition to achieve the DoF of KS. 

The following remark discusses the trivial case of the SVD-based OIA in terms of the antenna config- 
uration, where the inter-cell interference is perfectly canceled only through transmit beamforming. 

Remark 3: Note that if L > (K - 1)5 + 1, then G [iJl G c( K ~ 1)SxL in (HU) becomes a wide matrix 
and the singular value corresponding to v£ is 0. Therefore, w^l can be chosen such that rj^^ = 0. 
The result is intuitively immediate because the total rank of the effective interfering channels from each 
user to neighboring cells is (K — 1)5 and because at least one additional rank is required for each user 
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to transmit one data stream. From this result, it can be easily seen that in the case where each selected 
user transmits a < M data streams, all the inter-cell interference will be canceled through the SVD-based 
OIA if L > (K — l)Sa + a. In such case, the number of selected users, S, should be equal to or lower 
than [ M J. 



C. User Scaling Laws for Cell-Dependent L, M, N, and S 

In this subsection, we examine the user scaling laws for the case where L, M, N, and S are different 
from cells. Let us denote these parameters at the i-th cell by Li, Mi, Ni, and Si, respectively. The following 
theorem establishes the user scaling laws under this scenario. 

Theorem 3: With the cell-dependent parameters, the antenna selection- and SVD-based OIA schemes 
achieve KS DoF with high probability if 

Ni = u (SNR 5 ^) , and N, = u (SNR S ^" L ' +1 ) , (25) 

respectively, where S[ = £]JL fc=1 S k . 

Proof: See Appendix ITVl ■ 
From Theorem |3l it is seen that growing the number of serving users at the i-th cell, Si, increases the 
number of users required at all other cells for both the antenna selection- and SVD-based OIA. This is 
because increasing Si implies a reduced rank of the interference space at the i-th cell, on which users 
from the other cells attempt to align their signals. For the SVD-based OIA, large Lj reduces the user 
scaling condition of only the i-th cell. 



V. Comparison with Upper Bounds and Existing Schemes 

In this section, to verify the optimality of the proposed OIA schemes, we introduce an upper bound on 
the DoF. We also compare our schemes with existing schemes in terms of the achievable DoF and the 
computational complexity. 



A. Upper Bounds for DoF 

We now show an upper limit on the DoF in MIMO IMAC and discuss how to achieve the DoF upper 
bound. For completeness, we briefly review Corollary 1 of [21] in which the outer bound on the DoF of 
the MIMO IMAC is given by 

NK max (NL, (K - 1)M) 



DoF < min { NKL, KM, 



N + K-l 
NKmax((K- l)L,M) 
N + M-l 



(26) 



Now it is shown that choosing S = M, the proposed schemes achieve KM DoF with arbitrarily large 
N scaling according to (fT6l and (|24l) . Note again that with this choice, interference nulling is carried 
out through opportunistic user scheduling. As N increases, the outer bound (l2~6l is reduced to KM, and 
hence, our schemes can asymptotically achieve the optimal DoF. 



B. DoF Comparison with Existing Methods 

In this subsection, the proposed OIA schemes are compared with the two existing strategies [8], [21] 
that also achieve the optimal DoF in K-cell MIMO uplink networks. Let us first consider the K-uscr 
MIMO IC [8] with time-invariant or frequency-selective fading, which can be regarded as a MIMO IMAC 
with N — 1. Consider the case where L > M. Then both the scheme in [8] and the proposed SVD- 
based OIA with each user transmitting M spatial streams achieve the optimal DoF, given by KM, if 
L > KM [8, Theorem 3]. Note that in this case, interference can be perfectly nulled only through 
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SVD-based beamforming and thus no opportunistic gain is needed (refer to Remark [3]). The achievable 
scheme in [8] operates under time-varying or frequency selective fading channels with global CSI at 
all nodes, and the size of the time/frequency domain extension is given by (L/M + l)(n+ l) r , where 
T = KL/M ■ (K — L/M — 1) and n should be arbitrarily large to obtain KM DoF. For the if -user 
MIMO IC with time-invariant channel coefficients [8], a necessary condition for the parameter M is also 
needed to achieve the optimal DoF, which is given by M < (K — 2)L for K > 4. Hence, arbitrarily large 
M is also required as K increases, whereas our schemes have no necessary condition for M. 

Now, let us turn to the if -cell MIMO IMAC studied in [21]. For K = 2, both the transmit zero 
forcing scheme in [21] and the proposed SVD-based OIA with N = M achieves the 2M DoF if L > 
(K— 1)M + 1 = M + l. However, for K > 2, the scheme in [21] needs the necessary condition M > KS 
to obtain KS DoF [21, Theorem 3], which is not needed in the proposed OIA schemes. Moreover, the 
precoding matrices are designed based on the notion of global CSI in [21]. 

C. Computational Complexity 

In this subsection, we briefly discuss the computational complexity of the two proposed schemes and 
compare it to the complexity of the SISO IMAC scheme, Suh and Tse's scheme. The computational 
effort is analyzed in two-fold: the user computation and the BS computation. We omit the analysis of 
the detection and decoding complexity after the equalization at each BS, since it is all the same for the 
schemes considered. 

1 ) Antenna selection-based OIA: Each user calculates ([T3l , from which the scheduling metric (PT41) can 
also be obtained. From the results of [32], [33], it can be easily shown that the calculation of (fl"3l requires 
8(K — 1)MLS + Q(K — 1)LS — 2L floating point operations (flops), real additions or multiplications; 
thus, the complexity can be denoted by O(KLMS). 

Upon receiving N scheduling metrics, each BS selects S users with smaller scheduling metrics out of 
N users, which can be performed with linear-time complexity, i.e., 0(N), by the partial sorting algorithm 
[34]. Next, the construction of the effective channel matrix, i.e., F^ H (See ©), requires 8MS 2 — 2S 2 
flops. The inversion of this effective channel matrix to get F, needs 0(S 3 ) flops, and the calculation of 
r-j given in © requires 8MS + 8S 2 — AS flops. Therefore, noting that S < M, the overall computational 
complexity at each BS is 0{N + MS 2 ). 

2) SVD-based OIA: Each user first constructs defined in (fl"8~1) . which requires 8(K — 1)MLS — 
2{K — i)LS flops, i.e., O(KLMS). Note that the weight vector and the scheduling metric can be simulta- 
neously obtained from the SVD of G' M l The efficient and precise SVD method based on the Householder 
reflections and the QR decomposition can be performed with 0(KSL 2 ) flops [35]. Consequently, the 
computational complexity of the SVD-based OIA at each user is 0(KSL 2 + KLMS). 

All the procedure at each BS is the same as that of the antenna selection-based OIA except the construc- 
tion of the effective channel matrix, which requires 8MLS + 8MS 2 — 2MS — 2S 2 flops. Consequently, 
the overall complexity at each BS is given by 0(N + MLS + MS 2 ). 

Table I summarizes the computational complexity of the OIA schemes with the comparison to the 
SIMO case. It is obvious that the complexity is the lowest for the SEVIO OIA and is the highest for 
the SVD-based OIA. It is seen that as L increases, the complexity difference between the three schemes 
becomes greater. 

3) SISO IMAC: Now we briefly discuss the computational complexity of Suh and Tse's scheme [6]. 
Since this scheme applies only to the SISO IMAC, the comparison to this scheme is to roughly show 
the computational efficiency of the proposed schemes. Each user in Suh and Tse's scheme finds the 
inversions of K—l (nxn) -dimensional matrices and the Kronecker multiplications of K — 1 n-dimensional 
vectors, where n = K ~yN + 1. This calculation at each user requires 0(N + K K yiV ) flops. Another 
heavy calculation in this scheme is to find the (K — l)-level decompositions of n K ~ x x n A_1 matrices, 
which cannot be systematically performed. In addition, the complexity for the equalization at each BS 
is dominated by the effort to find the inversion of an {n K ~ l x n A_1 ) -dimensional matrix, which needs 
0(N 3 ) flops. 
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Considering the fact that Suh and Tse's scheme requires much lower dimension extension size than the 
conventional IA schemes and thus is already computationally attractive, the proposed schemes is more 
computationally effective compared to the previous schemes. In addition, it should be stressed that both 
the dimension extension size and N need to be arbitrarily large to achieve the optimal DoF with Suh 
and Tse's scheme, whilst arbitrarily large N suffices the condition for the optimal DoF for the proposed 
schemes. 

VI. Simulation Results 

In this section, through computer simulations, we evaluate the sum of LIF and the sum-rate of the 
proposed OIA schemes, operating with finite iV and SNR in the MIMO IMAC. The max-SNR scheme 
is compared, in which the weight vectors and the scheduling metrics are calculated at each user in a 
distributed manner only with local CSI. Specifically, each user employs eigen-beamforming to maximize 
its effective SNR and each BS selects 5 users having higher effective SNRs up to the 5-th largest one. 
The OIA scheme employing a fixed weight vector, i.e., = ei for all users, is also considered, which 
can be treated as the OIA scheme for SIMO IMAC. Thus, we refer this scheme as 'SIMO OIA. 

Figure ?? depicts the log-log plot of the sum of LIF, termed as sum-LIF, i.e., Ylf=i Sj=i versus 
N when K = 3, M = L = 2, and SNR is lOdB. This performance measurement enables us to measure 
the quality of the proposed OIA schemes, as shown in [5]. Specifically, Fig. ?? exhibits how rapidly 
the network becomes an error-free network with respect to N. Since the user selection of the max-SNR 
scheme does not contribute to the reduction of the LIF, the sum-LIF of the max-SNR scheme remains 
constant for increasing N. The sum-LIF of the antenna selection-based OIA decreases with respect to 
N at the same rate of the SIMO OIA, because the antenna selection-based OIA is subject to the user 
scaling condition SNR^ X_1 ^ 5 if L is fixed. On the other hand, the decreasing rate of the SVD-based OIA 
is higher, which is subject to the user scaling condition SNR ( ' X_1 ^ 5_L+1 . As 5 decreases, the decreasing 
rates of both the antenna selection- and SVD-based OIA schemes become higher due to the lowered 
scaling conditions. 

Figure ?? shows the log-log plot of the sum-LIF versus L when K — 3, M — 3, and N = 100. For 
the antenna selection-based OIA, the sum-LIF decreases linearly in log-log scale. On the other hand, the 
sum-LIF of the SVD-based OIA decreases much faster than the antenna selection-based OIA case and 
becomes zero if L > (K — 1)5'+ 1 = 5 (refer to Remark [3). Note that however, the feedback redundancy 
for the weight vectors grows as L increases in the SVD-based OIA, whereas no feedback is required 
regardless of L in the antenna selection-based OIA. 

Figure |6] depicts the sum-rates versus SNR when K = 3, and M = L = 2 for (a) N = 20 and (b) 
N = 100. The sum-rates of the considered schemes are saturated in the sufficiently high SNR regime, 
because the inter-cell interference cannot approach zero for fixed N values. That is, the SINR will be 
upper-bounded by a finite value for all schemes. In fact, 5 determines the amount of the interference 
level as well as the total DoF. For the max-SNR scheme, the interference at each BS increases as 5 
increases, whereas the sum-rate is increased by 5 times. The rate at each BS is approximately given by 
5 log ^1 + 1+ (j^?i)g.A ) ' where A denotes the amount of the interference received from a single user in 
a neighboring cell. Since this rate is a monotonically increasing function of 5, the rate of the max-SNR 
scheme grows with 5. On the other hand, the proposed schemes can significantly suppress the interference. 
Hence, the cases with 5 = 2 show higher sum-rates than the cases with 5 = 1 in the low SNR regime 
where the noise is dominant over the interference, and vice versa in the high SNR regime where it becomes 
more important to minimize the interference. As N increases, the interference can be more reduced, and 
thus the crossover SNR points, where the sum-rates for the cases 5 = 1 and 5 = 2 are identical, become 
higher. From Fig. |6] the crossover SNR points of the antenna selection-based OIA appear approximately 
at 6dB when iV = 20 and at 9.1dB when N = 100, whereas those of the SVD-based OIA are 8.1dB 
when N = 20 and 12.1dB when N = 100. 
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Figure [7J depicts the sum-rates versus N when K = 3, M = L = 2, and SNR is 20dB. For each of the 
scheme, the best S value was applied accordingly, which shows higher achievable rates. It is apparent that 
for infinitely large N, the rates of all the OIA schemes will be the same as those of the interference-free 
network. It can be seen from the figure that the SVD-based OIA with S = 1 approaches the upper-bound 
most rapidly, since the interference can be made smaller than that of the other OIA schemes according 
to the given scaling laws. While the SIMO OIA is inferior to the max-SNR scheme if iV < 20, both the 
proposed OIA schemes exhibit higher sum-rates than those of the max-SNR scheme if N > 3. 

Finally, Fig. [8] illustrates the symbol error rate (SER) averaged over all users versus iV when K = 3, 
M = L = 2, S = 1, and SNR is 20dB. The block length for each channel instance was assumed to 
be 50 symbols and quadrature phase shift keying (QPSK) modulation was used. For comparison, we 
considered the intercell interference-free scheme with the random user selection, which is labeled as 
'Interference-Free' in the figure. It is shown that the SERs of all the OIA schemes approach to the SER 
of the interference-free scheme as N increases. The trends for the approaching rates comply with the 
results of Theorem [T] and |2l that is, a lower user scaling condition implies better performance, a higher 
approaching rate in this case. 



VII. Conclusion 

We have proposed two OIA schemes for the MIMO IMAC and have derived the user scaling law 
required to achieve the target KS DoF. Although the antenna selection-based OIA cannot fundamentally 
change the user scaling law compared to the SIMO case, it can increase the achievable rate even with 
fixed L and with no feedback. Moreover, if L scales also with respect to SNR, then the scaling condition 
is linearly reduced with respect to L. It was also shown that the user scaling condition can be significantly 
reduced to SNR (iC-1) ~ i+1 using the SVD-based OIA with help of optimizing a beamforming vector at 
each user. Furthermore, the achievable rate of the proposed OIA techniques outperform the conventional 
user scheduling schemes including SIMO OIA. 

From this study on the user scaling law, we characterized the lower- and upper-bounds for the trade-off 
between the number of users required to achieve a target DoF and the amount of the feedback for the 
weight vectors. Even with the practical rages of the parameters, the user scaling law is a powerful tool 
to analytically compare the performance, such as the achievable rates or DoF, of any OIA schemes for 
given number of users. 

It can be conjectured that the MIMO OIA with limited feedback for the weight vectors will make a 
bridge between the proposed two OIA schemes. As our future work, the scaling law for the number of 
users as well as the feedback size will be studied. 



Appendix I 
Proof of Theorem Q] 

From CEO) and (HI), SINR^' 1 can be written as 

SNR 



SINR 



■*j i 



*i,3 



> 



SNR/ ||t 



1,3 



,k^i Z—dm=l 



(27) 
(28) 



i,l(k,m) 



SNR 
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It is apparent that the DoF of KS is achieved if the interference term in the denominator of the right-hand 
side of (|28l) remains constant for increasing SNR. At this point, let us define Pas by 



V, 



AS 



K S 

lim Pr<( V Vllu^h 1 ' 

k=l,k^i m=l 



SNR->oo 



Hy,[k,m] 

l(k,m) 



SNR < e, 



V user j in the i-th cell, i E KZ,j G S 



(29) 



where e > is a positive constant. Then, DoF is bounded as 

DoF > KS ■ Pas- 



(30) 



When calculating the lower bound (|30l ), we assumed that the DoF of KS is achieved if the interference 
remains constant for increasing SNR, and zero DoF is achieved otherwise. 

The essential of the OIA is the fact that the sum of the received interference terms is equivalent to the 
sum of the LIF metrics of the selected users. That is, 



i,l(k,m) 



EE*" 

i=i j'=i 



K 



H u [k,m] 

i,l(k,m) 



K K S K S 

E E Ellu.^ 

i=l fc=l,t/i m=l 

Subsequently, defining 

E Ell u * h 

k=l,k^i to=1 

we find the following lower-bound of Vas- 

{K S 
VV/ ASii SNR<e 
i=l j=l 

{S K S 
VV Yva( ] SNR<e 
j=l i=i j>=i 

> lim Prl^f < SI ! R ~ 9 l£ ,Vi e /c,v/ e s 



(31) 



(32) 



SNR^oo 



KS 2 



(33) 

(34) 
(35) 



' AS 



where (|34l) follows from (I31D . Unlike in the SIMO case [16, Theorem 1], rj A i! is the minimum of L 



H u [kM 



, I 



independent Chi-square random variables with degrees-of- freedom of 2(K — 1)S, 
1, . . . , L. We denote the probability that user j in the i-th cell has at least one transmit antenna with the 
scheduling metric lower than £S ^ 2 as 



^ = i-H E || U ^ h M 

k=l,k^i 



2 eSNR 

> 



KS 2 



V/G{1,...,L} 



(36) 



It can be easily verified that P a is identical and independent for all users. Let us denote the right-hand 
side of (1331 ) by Vl s . Note that V AS represents the probability that there exist at least S users in each cell, 
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which have the scheduling metrics lower than 6 ^ 2 , and thus we have 

s-i / M \ 

v° AS = i- lim y; ( 7 ) p: . (i - p a f-\ ov) 

SNR^oo-^^ V I I 

i=0 V 7 

Denoting by F(x) the cumulative density function (CDF) of a chi-square random variable with the degrees- 
of-freedom of 2(K — 1)5, we have 

/ r.MNK ' \ \ L 

P a = l-[1-F 



Applying (|38l ) to 071 ), we get (|39l and (|40l ) at the bottom of the next page, where C\ and C 2 are constants 
independent of SNR and L, defined by 

e -i2~(K-i)s , £ wi<:~i)s 

Cl = (K -1)S -T((K -1)S) ' \~KS*) ' (41) 

2 -(*-i)s+i / e 

^ = (K-l)S-T((K-l)S) ' \YS^) ' (42) 



< F ^ (43) 



Here, (|40l) follows from the fact that [15, Lemma 1] 

e _l 2 _(^_l) 5 



- i)^ • r((K - 



o-(if-l)S+l 

Ffx) < - x' {K - 1)s (44) 

1 ' ~ (K-1)S-T((K-1)S) { ' 

and from the fact that ^^l^, < N l . Here, if we choose e small enough such that C 2 SNir (A '- 1)5 < l/L 
for given SNR, we get 

(l - C 2 SNR- {K - 1)s ^j > 1 - LC 2 SNR- {K ~ 1)S , (45) 

which follows from the fact that 1 — xy < (1 — x) v for any < x < 1 < y and xy < 1. Now, inserting 
(1431) to (001) gives us 

P AS >1- lim V-i 1 1 — (46) 

as- snr-,^ (1-C 2 SNR 5 ) LJ 

where 5 = -{K - 1)5. If LN = co (SNR {A '~ 1)S ) , then (l - dSNR"^ 1 ^)^ decreases exponentially 

with respect to SNR, whereas ^iVLC 2 SNR~ ( - A ~ 1 ' )S j increases polynomially for any i > 0. Therefore, 

Vl s tends to 1 as SNR goes to infinity, and thereby Vas tends to 1. This proves the theorem together 
with d30h. 



P° s = 1 - lim V — '- n (39) 



s-i (n ( 1 - ( 1 - C 2 • SNR-^- 1 ) 5 ) U f 1 - dSNR-t*- 1 ^ 

>1- lim V 1 V V )l \ — , (40) 

snr— x / -— ' (l - C 2 SNR-( A -^ 
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Appendix II 
Proof of Lemma Q] 

H 7*1 

Since is chosen from an independent isotropic distribution and rl k is an i.i.d. complex Gaussian 
random matrix, for all i, k E /C, j E S, is also an i.i.d. complex Gaussian random matrix. 



Furthermore, both of and Hj! are chosen from the continuous distributions, and thus have full 
ranks almost surely [36]. The LIF metric = oj*'^ is the smallest eigen value of the (L x L)- 

dimensional central Wishart matrix G^ G^. Therefore, from [37, Theorem 4], the polynomial CDF 
of the smallest eigen value of the full-rank Wishart matrix which is constructed from a ((K — 1)S x L)- 
dimensional complex Gaussian matrix has the smallest power of (K — 1)S — L + l with the multiplicative 
coefficient a defined by 

rwi) 



a: 



((K-1)S-L + 1)\T L (L) 1 1 

Here, T s (t) is the normalized complex multivariate gamma function, i.e., T s (t) = Yli=i{t ~ 
an (L x L) -dimensional integer matrix defined as 

L — i\ i — 1, 



'■j 



.,L- l,j = 1, . . .,L, 

j > i 

,L,j = l,...,L,j < i 
otherwise. 

Therefore, a is determined only by K, S, and L, which proves the lemma. 



J - « 

(-iy-HL-j)\ 
(n-jy. 





Appendix III 
Proof of Theorem [2] 
From the SINPJ M ' lower bound, given by 

SNR/llf^H 2 



SINR l '" iJ > 



1 + Ylk=l,k^i Hm=i 

we again consider the lower bound of the DoF as 

DoF > KS ■ Psvd 



U « 1 1 W SVD 



SNR 



(47) 
and S is 

(48) 



(49) 



(50) 



P S vd= lim Pr^ IsmiSNR < e, 

SNR^oo 



V user j in the i-ih. cell, i E IC,j E S 



(51) 



A' 



flTjl'i" 1 ] [k,m] 



W 



SVD 



where 

^svd,* = ^2 || u * Hh !' 

k=l,k^i m=l 

Similarly to (|3TT > to (1331 ), the lower bound on Psvd is obtained from 

{K S 
VV/svD^SNR^e 
»=i i=i 



I j=l i=l j'=l J 
> lim p r (4«5<^^!,Vz6X:,Vi / e<S 



SNR^oo 



(52) 

(53) 

(54) 
(55) 
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The right-hand side of (1551 ) is the probability that there exist at least S users with the scheduling metrics 
lower than SN ^ S i e ■ Noting that the scheduling metrics 77svr)' i — ^, ■ ■ ■ , K, j — 1, . . . , S, axe identically 
distributed, the right-hand side of (1551 ), denoted by V$ YD , can be expressed as 



i=0 



Denoting p = (K - 1)S + L — 1, we further have 

51 iV! 
"po = 1 - lim > — — - 

SVD SNR^oo^ i\(N-i)\ 



(^SNR^ + o(SNR^)) J 
(1 - ^SNR" P - o (SNR" p )) i 
x (1 - ^SNR~ P - o (SNR-"))* (57) 

> 1 - lim > i — i i '-^-r- 

snr^oo^ (l -$SNR- p -o(SNR- p ))* 

x (1 - ^SNR~ P - o (SNR _P )) 7V (58) 



where 



/ e s(K-l)S-L+l 



Here, (|57l) follows from Lemma [T| and from choosing e small enough such that eS ^ 2 1 < 1 for given 



SNR, and W follows from < N l 

Now, if V = to (SNR P ), (l - ^SNR~ P - o (SNR~ P ) ) N decreases exponentially as SNR increases. On 
the other hand, {N (\I/SNR _p + o (SNR~ P )) }* increases polynomially for any i > 0, and thus, the second 
term of (1581) tends to zero as SNR — ► oo. Therefore, the lower bound of "Psvd given in (1551) tends to 1, 
which proves the theorem together with (|50T ). 

Appendix IV 
Proof of Theorem [3] 

Following (0 to d35j) and <|53]) to (O and denoting the scheme indicator by r G {AS, SVD}, V T with 
the cell-dependent parameters can be written by 

^ > s N feoo Pr IE E E^" ] SNR < e | (60) 

I i=l i=i j'=l J 
SNR x e 

SNR— >oo I - Ef=l^' 



> Jim Pr <| i 7 W'] < ^ — ,Vi G /C, V/ G S \ (61) 



— S^r'S'"^^' <62) 
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where in (1621 ), Vr denote the probability there exist at least Si users with LIF metrics smaller than 



snr- 



at the i-th cell, which is independent from those of the other cells. 



i) Antenna selection-based OIA: Since 



U fc "h 



k,l(i,j) 



is a Chi- square random variable with DoF 



of 2Sk, the scheduling metric 77^5 in ([141) is a Chi-square random variable with DoF of 25", where 
S' = 2 Ylk^i k=i The rest of the proof can be done analogously to the proof for Theorem [TJ replacing 
(K - 1)S with S'. 

ii) SVD-based OIA: Since U k H W' j] is an (Sfc x Li) dimensional Gaussian matrix, defined in 

(fl"8l) is now (S' x L,)-dimensional. Following the analogous derivation of the proof for Theorem |2] and 
replacing (K — 1)S with S', we can complete the proof. 
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TABLE I 

Computational complexity of the OIA schemes (flops). 
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Fig. 3. Proposed MIMO OIA where K = 2, M = 3, and S = 2. 
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Fig. 4. Sum-LIF versus N for the MIMO IMAC with K = 3 and M = L = 2. 
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Fig. 5. Sum-LIF versus L for the MIMO IMAC with K = 3, M = 2, and N = 100. 
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Fig. 6. Achievable sum-rates versus SNR when K = 3, M = L = 2, and (a) N = 20 and (b) TV = 100. 
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Fig. 7. Achievable sum-rates versus TV when K — 3, M = L — 2, and SNR=20dB. 
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Fig. 8. Average SER versus N with QPSK signaling when K = 3, M = L = 2, S = 1, and SNR=20dB. 



